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FOREWORD 


The  research  work  described  in  this  report  was  performed  by  the  Bell  Aircraft  Corporation, 
Buffalo,  New  York,  for  the  Aircraft  Laboratory,  Wright  Air  Development  Center,  Wright-Patterson 
Air  Force  Base,  Ohio,  The  work  was  accomplished  under  Air  Force  Contract  No.  AF33(616)-2536, 
Project  No.  1370,  (UNCLASSIFIED  TITLE)  "Aeroelasticity,  Vibration  and  Noise,"  and  Task  No. 
13479,  (UNCLASSIFIED  TITLE)  'Investigation  of  Flutter  Characteristics  of  All  Movable  Tails." 
Henry  Katz,  1st  Lt,  USAF,  of  the  Dynamics  Branch,  Aircraft  Laboratory  is  task  engineer.  The 
transonic  program  was  carried  out  by  the  Aeroelasticity  Section  of  the  Bell  Aircraft  Corporation 
Research  Division,  with  Mr.  B.  B.  D' Ewart,  Jr.  acting  as  project  engineer,  assisted  by  Mr.  R.  F. 
Farrell.  Research  started  in  May  1954  and  is  part  of  a  continuing  program  of  aircraft  flutter 
research  by  WADC. 

This  document,  excepting  the  title,  is  classified  CONFIDENTIAL  because  lt  reveals  information 
which  can  be  employed  to  determine  design  criteria  for  the  prevention  of  flutter  of  all- movable  tails 
at  trar^oaic  speeds. 
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ABSTRACT 

\ 


Wind  tunnel  tests  were  conducted  on  a  family  of  10  all  movable  45-degree  swept  stabilizer 
models.  The  model  panels  were  essentially  dynamically  Identical  except  that  panel  stiffness  levels 
Were  varied  to  obtain  flutter  within  the  Mach  0.8  to  1.2  speed  range  at  a  mass  ratio  of  25.  For  pitch 
to  bending  frequency  ratios  of  2.67  and  1.50,  and  a  mass  ratio  of  25,  flutter  speeds  remained  between 
16  and  19  percent  higher  than  predicted  by  incompressible  swept  analyses  using  the  methods  of  NACA 
TR1014.  No  significant  transonic  dip  or  supersonic  rise  in  experimental  values  of  VE/br  uja  was  ob¬ 
served  for  the  above  conditions,  however,  a  region  of  low  stability  was  found  to  exist  below  the  flutter 
boundary  beyond  Mach  1.0. 

The  effects  of  variations  in  pitch  to  bending  frequency  ratio  with  mass  ratio  constant  at  25.0 
are  small  and  of  the  same  order  of  magnitude  as  *re  predicted  by  reference  analysis. j^pata  points 
below  Ma£lr~4h9^_at  u  =  25,  show  a  drop  of  5%  in  VE/br  tuhj  as  UJ^/aihj  ^drbps  from  2.67 
to  0.9  Analyses  shovTa~co«£gpondlng  drop  of  8%.  The  drop  in  V'E/br  UJh,  in^getngfrom  co^/cuh,  = 
2.67  to  1.50  does  not  exceed  ^bTirany-pQiat_between  Mach  0.8  and  1 .  ljj.at"/x=  25.0.  Increasing  mass 
ratio  does  not  bring  about  as  much  of  anTryrgrmr  in  by  incompressible 

analysis  particularly  at  cu^/cbh,  =  0.9.  At  this  frfiguenCyratio  at^icTri>19  and  /x  =  70  the  in¬ 
compressible  analyses  gave  a  reduced  velocify-tilSfwas  10  percent  unconservativfe.  ^Test  data  were 
not  sufficient  to  determine  the  degree  of  interdependence  of  Mach  number  and  mass  ratio  effects. 

At  Mach  0.8  the  flutter  to  torsion  frequency  ratio  Oig/wa  was  from  0.55  to  0.60  and  was  little 
affectedby  variations  in  cu^/u^,  .  Increasing  Mach  number  up  to  1.2  at  yx  =  25  resulted  in  a  drop  in 
CJE/o;a  of  about  20  percent.  Increasing /x  above  25  causes  further  drops  in  G0E/aja  .  For  the 
configuration  tested  pitch  does  not  become  significant  amplitude-wise  until  cox/cuf,  drops  down  close 
to  unity.  Y 
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INTRODUCTION 


There  has  been  a  definite  trend  towards  utilization  of  all-movable  control  surfaces  on  modern 
high  speed  aircraft  to  obtain  increased  control  effectiveness.  However,  the  use  of  this  type  of  control 
clearly  increases  the  possibility  of  flutter  problems  as  is  substantiated  by  the  fact  that  at  least 
seven  cases  of  flutter  of  this  type  of  surface  are  known  to  have  occurred  on  U.S.  aircraft  and 
missiles  within  the  last  four  years.  This  has  pointed  up  the  need  for  a  better  understanding  of,  and 
control  over,  the  parameters  contributing  to  flutter  of  these  surfaces.  This  report  summarizes 
wind  tunnel  test  data  in  the  Mach  0.8  to  1.2  speed  range  on  a  family  of  10  all  movable  stabilizer 
models  of  45  degree  quarter  chord  sweep.  Model  geometry,  mass  and  elastic  characteristics  are 
summarised  in  Section  1  of  this  report.  These  characteristics  are  typical  of  existing  aircraft  con¬ 
figurations.  Design  and  construction  details  are  presented  in  Appendix  L  Model  calibration  and 
data  reduction  procedures  are  covered  in  Appendix  II. 

Testing  was  accomplished  using  the  WADC  10ft.  Transonic  Tunnel.  Models  were  mounted  from 
a  side  wall  of  the  test  section. 

Flutter  analyses  were  performed  using  incompressible  two-dimenBional  air  force  coefficients 
in  analyses  utilizing  the  methods  of  NACA  TR  1014  (Ref.  1).  Analysis  results  are  summarized  in 
Section  2.  Presentation  of  flutter  test  results  is  made  in  Section  3,  A  discussion  of  these  results 
and  a  comparison  of  analysis  and  test  results  is  also  made  in  Section  3. 


Manuscript  released  by  the  Authors  Apri  I 
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SECTION  I 

MODEL  CHARACTERISTICS 


The  general  features  of  the  test  configurations  are  shown  In  Figure  1.  The  10  models  tested 
were  geometrically  identical  and  had  similar  mass  and  elastic  characteristics  except  that  model 
stiffness  levels  were  varied  to  allow  flutter  within  the  Mach  0.8  to  1.2  speed  range  at  a  mass  ratio 
of  25.  Flexure  pivots  provided  a  pitch  degree  of  freedom  through  the  35.0%  MAC  based  on  span  to 
fuselage  centerline.  By  using  two  flat  torsion  springs  with  variable  effective  lengths,  pitch  frequency 
could  be  varied  from  33  cps  to  214  cps.  Data  was  taken  almost  entirely  at  pitch  to  bending  frequency 
ratios  of  0.9,  1.5  and  2.87.  A  method  was  provided  for  locking  out  the  pitch  degree  of  freedom  to 
allow  the  stabilizer  models  to  be  checked  as  cantilevers.  General  model  characterir  ti.cs  are  summar¬ 
ized  beneath  Figure  1.  Important  parameters  of  the  individual  models  are  summarized  in  Table  I. 
Model  dimensions,  construction  details,  etc.  are  covered  in  Appendix  1. 

'  A  half  fuselage  of  5-inch  radius  was  provided  to  reduce  tunnel  boundary  layer  effects.  The 
test  assembly  was  mounted  from  a  side  wall  of  the  WADC  10  ft  Transonic  Tunnel.  One  of  the 
tunnel's  circular  windows  was  replaced  with  a  steel  plate  to  which  the  model  supporting  structure, 
operating  mechanisms,  and  fuselage  section  were  attached.  The  window  frame  contained  a  drive 
mechanism  which  allowed  the  rotation  of  the  entire  assembly  for  streamlining  the  models. 

Models  were  provided  with  bending-sensitive  and  torsion-sensitive  strain  gage  bridges  for  use 
in  determining  flutter  frequencies  and  amplitude  ratios.  The  flexure  pivots  in  the  pitch  system  were 
also  strain  gaged  to  provide  corresponding  data  for  pitch.  Details  of  model  construction  and  cali¬ 
bration  are  given  in  Appendix  n. 


WADC  TR  57-392  2 

UNCLASSIFIED 


UNCLASSIFIED 


Figure  1.  General  Model  Totting  Configuration  Used  in  WAOC  10  Ft  Transonic  Tunnel 


MODEL  CHARACTERISTICS; 

45.  Deg, 

0.5 
2.91 

NACA  65  A 004  In  Wind  Direction 
35.0%  MAC  Based  on  Span  to  Fuselage 
Centerline 

or 

60.4%  Root  Chord  in  Wind  Direction  at 

Surface  of  Fuselage 

E.A.  41.6%  Chord  Normal  to  Quarter  Chord 
c.g.  50%  Chord  Normal  to  E.A, 

For  detailed  dimensions  see  Appendix  I 

Range  of  Torsion  Frequencies  of  10  Models  Tested,  153  to  219  cpe 
Average  Sending  To  Torsion  Frequency  Ratio  0.273 
Variable  Pitch  Frequency  Range  33  to  214  cps 
For  Masa  and  Elastic  Characteristics  of  Each  Model  See  Table  I 


Quarter  Chord  Sweep 
Taper  Ratio  of  External  Span 
Aspect  Ratio  of  External  Span 
Airfoil 

Pitch  Axis  location: 
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Averages _ _ _  .273  20.53  .168  .248  41.6  5 

Pitching  Inertia  of  all  elements  of  pitch  system  behind  fuselage  fairing  «  2.05  lb  in^  (varies  slightly  with  effective  pitch  spring 
length)  see  Fig.  25 
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SECTION  2 

FLUTTER  ANALYSES 


Four  degree  of  freedom  flutter  analyses  were  performed  to  provide  the  usual  incompressible 
flow  "reference  analysis"  data  for  use  as  a  basis  of  comparison  of  analysis  versus  experiment.  First 
bending,  torsion,  second  bending  and  pitch  degrees  of  freedom  were  used.  Uncoupled  pitch  to  bending 
frequency  ratios  and  mass  ratio  were  varied  over  the  range  of  conditions  encountered  in  the  test 

program. 

Analyses  were  performed  using  for  xa  the  original  design  value  of  0.20.  Since  laboratory 
tests  showed  variations  in  xa  down  to  0.15  a  number  of  check  points  for  xa  =  0.17  and  0.15  were 
made.  Resulting  variations  in  VR/br oja  were  1.0%  or  less. 

Development  of  the  flutter  equations  follows  the  general  methods  outlined  in  NACA  Report 
1C14  {Reference  1)  expanded  to  include  the  pitch  and  second  bending  modes.  Secondary  spanwise 
tiow  terms  as  indicated  in  Reference  1  by  heavy  brackets  have  been  omitted  in  accordance  with 
common  practice.  Analysis  details  and  a  summary  of  determinant  elements  are  given  in  Appendix  m. 

Analysis  results  are  summarized  in  Figures  2  through  6  which  show  the  flutter  parameters  as 
functions  of  mass  ratio  and  cu£/ajh|  .  Figure  2  shows  the  variations  in  Vp/ty  gl>0  with  mass  ratio  and  pitch 
to  bending  frequency  ratio.  Goingfrom  /x=  25  to  fi-  75  at  constant  cux/Wh,  causes  a  43  to  54  percent 
Increase  InVfa/ty&bdepending  on  the  value  of  cu^/cuh,  while  a  change  in  gux/oj^  from  0.9  to  oo  at  constant 
causes  only  a  10  to  17  percent  increase  in'^/brU)a  depending  on  the  value  of  u. .  Figure  3  shows 
that  a  change  in  mass  ratio  from  25  to  70  while  holding  Oty/coh,  constant  causes  a  16  to  26  percent 
drop  In  <uR/oua  depending  on  the  value  of  oj^/cuh|  ,  while  a  change  in  c from  0.9  torn  at  constant 
fi  causes  a  26  to  44  percent  increase  in  OJp/ cl) a  . 

The  reference  analyses  show  flutter  for  a  given  model  ( tyuh(  =  a  constant)  as  a  f  -notion  of 
true  airspeed  V  and  mass  ratio  .  Selecting  V/tv  ^x*  as  a  flutter  index,  x  may  be  solved  for 
to  make  this  Index  a  constant  between  specified  values  of  p  .  Experience  to  date  indicates  that 
x  will  lie  between  0.5  and  0.0.  When  x  =  0.5,  the  index  becomes  the  familiar  V/tvcdh,  which  has 
the  special  feature  of  being  a  function  of  q  .  Thus,  a  constant  value  of  this  flutter  index  for  a  given 
model  or  airplane  indicates  a  constant  q  flutter,  x  approaching  zero  indicates  decreasing  sensitivity 
of  flutter  speed  to  changes  in  mass  ratio,  with  velocity  becoming  invariant  at  x  =  0  for  a  given 
model.  Since  stability  then  depends  only  on  true  airspeed,  this  limiting  condition  has  been  called 
a  true  airspeed  flutter.  For  a  family  of  models  differing  only  in  natural  frequencies,  V/bra>h  rather 
than  V  would  be  Invariant  for  the  true  airspeed  flutter  condition. 

Solving  for  x  as  indicated  above  using  the  data  of  Figure  2  over  increments  of  mass  ratio  of 
10.0  gave  values  of  x  ranging  from  0.47  at  (J.  =  25  range  down  to  0.28  at  ft  =  70,  for  ai  x  /ojh|  =  co 
and  corresponding  values  of  x  irom  0.42  to  0.20  for  u)^,/aih|  =  0.9.  Figure  4  shows  the  variations 
in  VR/br  tUh,  V/x over  the  /x  range  used  in  testing.  If  analyses  predicted  corstar.t  q  flutter,  the 
curves  of  Figure  4  would  be  straight  horizontal  lines.  A  constant  q  flutter  would  show  up  on  Figure 
2  as  Indicated  by  the  dotted  line  for  the  =  co  case  based  on  the  value  of  VR/br  aih|  at 

/x  =  25.  Figure  5  shows  the  drops  in  q  at  flutter  relative  to  the  value  at  fx  -  25  which  occur  as 
/x  is  varied  from  25  to  70. 

To  show  the  small  sensitivity  of  VR/br  uh|  to  changes  in  pitch  to  bending  frequency  ratio, 
Figure  4  has  been  crossplotted  to  give  Figure  6. 
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Figure  2.  Analysts  results,  Vr  /br  u>a  versus  mass  ratio  for  pitch  to  bending  frequency  ratios 
used  in  testing.  Dotted  line  indicates  constant  q  flutter  (Le.VR/br<oa «  a  (constant)  with 
=  ®  >  based  onVR/b,  ciij,  ot/i.  ■  25 
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bending  frequency  ratios  used  in  testing. 
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Figure  5.  Analysis  results,  variations  in  q  as  a  ratio  with  the  value  at  /a  =  25  versus  mass 
ratio  for  pitch  to  bending  frequency  ratios  used  in  testing. 
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SECTION  3 
TEST  RESULTS 


1.  Presentation  of  Data 

A  summary  of  test  data  Is  given  In  Table  IV.  With  one  exception  test  points  were  limited  to  pitch 
to  bending  frequency  ratios  of  0*9.  1.5, 2.67  and  ®  .  The  data  at  these  frequency  ratios  is 

plotted  on  Figures  7  through  10  in  terms  of  the  parameter  VE/b,  o»h,  JJx  .  As  much  data  as  possible  was 
obtained  at  a  mass  ratio  of  25  to  eliminate  this  parameter  as  a  test  variable.  Tests  were  conducted 
at  high  mass  ratios  In  some  cases  to  get  flutter  at  high  Mach  numbers  where  frequencies  of  available 
models  were  too  low  to  give  the  desired  test  speed  at  fx=  25.  In  such  cases  the  desired  Mach  number 
was  held  and  test  density  was  increased  to  get  flutter.  Figure  11  shows  all  data  in  terms  of  the  re¬ 
duced  velocity  parameters  Ve/bfO)h|  and  Vg/b,  OJa  i 

In  considering  the  high  mass  ratio  points  on  Figure  11  it  should  be  realized  that  at  a  given  Mach 
No.,  the  ordinates  are  a  function  of  ,  ft  ,  and  uty/cuh,  •  Therefore  conclusions  as  to  the  effects 
of  any  one  of  toe  variables  can  be  made  only  for  those  limited  cases  where  the  other  variables  were 
held  constant.  Ratios  of  experimental  flutter  velocities  to  corresponding  analysis  velocities  were 
obtained  using  the  analysis  data  of  Section  2.  These  velocity  ratio  data  are  shown  graphically  in  Fi¬ 
gure  12.  The  average  value  of  the  speed  of  sound  at  the  test  section  was  1060  ft/sec  and  the  average 
deviation  from  this  value  was  1.2%. 

A  high  degree  of  tunnel  turbulence  plus  the  large  region  of  very  low  stability  below  the  true 
flutter  boundary  in  the  transonic  region  made  the  identification  of  the  true  flutter  boundary  very 
difficult,  time  consuming  and  costly  in  models,  and  was  a  major  factor  in  limiting  the  total  number 
of  data  points  obtained.  Sufficient  data  were  obtained  at  pitch  to  bending  frequency  ratios  of  1.5  and 
2.67  at  a  mass  ratio  of  25  to  allow  trend  curves  to  be  drawn  for  the  respective  V^/br  <Oh,  SpL  plots  of 
Figures  8  and  9.  /x=  25  data  at  these  two  frequency  ratios  on  Figure  10  through  12  has  been  repre¬ 

sented  by  the  trend  curves  rather  than  by  the  data  points.  Data  at  other  frequency  ratios  occurred 
over  such  a  range  of  mass  ratios  and  were  so  sparse  that  trend  curves  were  not  attempted.  The 
separate  data  points  are  therefore  shown  and  individually  identified  on  Figures  10  through  12. 

Flutter  frequencies  are  shown  ratioed  with  their  respective  torsion  frequencies  and  analyses 
flutter  frequencies  in  Figures  15  and  16  respectively.  Amplitude  ratio  data  are  presented  in  Figure 
17. 


2.  Discussion 

a.  Dat4  at  Oty/wh  =  0.9 

Since  date  at  this  frequency  ratio  were  considered  of  less  practical  importance  for  present 
day  aircraft  Chan  date  at  higher  frequency  ratios,  tests  at  this  frequency  ratio  were  given  a  low 
priority.  As  a  result  only  3  flutter  points  were  obtained  as  shown  on  Figure  7.  The  large  increases 
in  the  mass  ratios  (necessary  to  get  models  8  and  10  to  flutter  supersonically  at  this  low  pitch  to 
bending  frequency  ratio)  caused  a  pronounced  drop  in  vE/br  cuh,  as  the  Mach  number  was 
Increased. 


Test  points  at  Macli  0.80  and  1.05  involved  a  simultaneous  change  in  3  variables:  cun,  , 
Mach  number  and  fM  ,  making  trend  observations  difficult.  Th  'ach  1.05  and  1.19  test  points  were 
obtained  using  the  same  mode.,  and  therefore  these  two  test  j  may.be  considered  as  showing  the 
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effect  of  a  change  in  speed  and  altitude  on  a  given  stabilizer.  Going  from  Mach  1.05  to  1.19  repre¬ 
sents  a  13.3%  increase  in  Mach  number.  However,  due  to  a  33 °F  drop  in  tunnel  stagnation  tempera¬ 
ture  plus  the  increased  temperature  drop  due  to  the  expansion  to  the  higher  Mach  number,  the  cor- 
respon  ling  increase  in  velocity  and  hence  V^/br  aih,  was  only  2.0%.  This  indicated  velocity  increase 
combined  with  the  change  in  mass  ratio  from  40.1  to  70.9  resulted  in  a  40%  drop  in  dynamic  pressure. 
The  reference  analysis  predicts  a  17%  increase  in  Vp/br  uih|  f°r  315  increase  in  jj,  from  40.1  to  70.9%. 
The  failure  of  this  large  increase  in  reduced  velocity  to  show  up  experimentally  points  toward  an 
unconservative  tendency  in  "reference"  analyses  as  [±  is  increased  and/or  the  existance  of  a  large 
destabilizing  compressibility  effect  acting  at  this  pitch  to  bending  frequency  ratio.  Additional  data 
points  at  [±  =  25,  40  and  70  throughout  the  entire  high  subsonic  and  transonic  speed  range  wouia  oe 
valuable  in  isolating  compressibility  and  mass  ratio  effects. 

b.  Data  at  cu^/Wh,  =15 

At  a  mass  ratio  of  25  a  flutter  boundary  has  been  established  as  shown  on  Figure  8.  For 
this  mass  ratio  and  frequency  ratio  the  possibility  of  a  "transonic  dip"  between  Mach  0.95  and  1.15 
(where  actual  flutter  points  are  lacking)  is  eliminated  by  virtue  of  the  test  sequence  used  in  arriving 
at  the  Mach  1.18  test  point  (model  6,  f±  =  28.0).  The  test  path  is  plotted  on  Figure  8  as  a  phantom 
line.  It  is  seen  to  approach  the  [±=25  flutter  boundary  along  a  tangential  path  through  the  transonic 
region.  Since  the  mass  ratio  along  this  path  was  always  in  the  vicinity  of  25,  any  pronounced  "tran 
sonic  dip"  would  have 'caused  flutter  to  appear  at  a  considerably  lower  Mach  number  than  was  found 
to  be  the  case.  Reference  3  reports  a  small  dip  for  similar  swept  cantilever  models  at  [±  =25.  The 
data  point  obtained  on  model  6  at  Mach  1.18  was  a  well  defined  divergent  flutter.  Figure  13  is  a 
section  of  oscillograph  record  from  this  run  just  before  the  model  was  lost.  Since  it  occurred 
at  virtually  the  same  value  of  Vg/br  oih,  </f±  as  did  the  subsonic  data  points,  no  significant  super¬ 
sonic  rise  of  the  type  reported  in  Reference  3  for  similar  swept  cantilevered  models  seems  to  exist 
at  [±=25  for  this  pitch  to  bending  frequency  ratio  and  Mach  range. 

A  very  low  stability  region  below  the  flutter  boundary  was  found  to  exist  above  Mach  1 .0 
for  [±  greater  than  25.  This  has  been  Indicated  in  Figure  8  as  a  dotted  area.  Models  in  this  region, 
always  appeared  to  be  very  close  to  flutter  while  mass  ratios  were  always  above  25  and  sometimes 
above  50.  Even  though  models  were  near  flutter  throughout  this  region  it  was  usually  necessary  to 
make  large  changes  in  mass  ratio  to  cause  actual  flutter,  thus  indicating  a  low  sensitivity  to  density 
clianges. 

Comparing  Figure  8  with  Figure  7,  it  appears  that  the  low  stability  at  oj^/cuh,  =1.5  shown  in 
Figure  8  is  altered  to  a  condition  of  instability  when  u>^/u;hj  is  dropped  to  0.9  causing  flutter  to 
occur  along  the  lower  side  of  what  was  the  low  stability  regions  for  cui/uih,  =  1.5.  The  [±=  45.4  data 
point  obtained  on  model  7  at  Mach  1.05  was  obtam^d  by  holding  Mach  number  constant  and  slowly 
increasing  test  density  (moving  straight  up  on  Figure  8).  The  model  fluttered  to  destruction  at  the 
point  indicated.  The  [±  =  38.9  data  point  shown  for  model  10  at  Mach  1.18  was  obtained  in  the  same 
manner;  however,  the  test  was  stopped  before  the  model  attained  a  divergent  flutter.  Going  from  this 
data  point  to  the  Model  6  data  point  at  Mach  1.18  and  [±  =  28.0  represents  a  drop  in  [±  of  2G%  with  an 
accompanying  drop  ir.  Vg/brouh,  of  only  4%.  A  minor  drop  in  reduced  velocity  (increased  stability 
due  to  increased  u)(,  )  is  equivalent  here  to  a  large  drop  in  mass  ratio  (decreased  stability  due  to 
larger  aerodynamic  forces).  Solving  for  the  value  of  x  which’ makes  the  flutter  index  \%/br  cuhi  [±l 
a  constant  for  these  two  flutter  points  gives  x  =  0.11.  The  two  flutter  points  therefore  indicate 
flutter  in  this  test  region  to  be  more  nearly  a  function  of  true  airspeed  than  q,  since  constant  q 
flutter  requires  that  x  =  0.5  in  the  above  flutter  index  (see  bottom  half  of  Page  5).  These  two 
test  points  and  other  high  mass  ratio  points  are  also  shown  on  the  Vjt/Vr  plot  of  Figure  12  illus¬ 
trating  the  decreas  in  conservatism  of  "reference"  analysis  relative  to  test  as  mass  ratio  is  in¬ 
creased.  It  appears  that  this  loss  of  conservatism  of  "reference"  analyses  with  increasing  mass 
ratios  is  due  to  the  relatively  small  increase  inN/e/brOth,  observed  experimentally  with  increasing  [± 
as  compared  with  reference  analysis  predictions.  While  experimental  data  is  insufficient  to  prove 
that  contributions  from  compressibility  and  pitch  to  bending  frequency  ratio  are  not  significant  In  this 
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connection,  such  a  conclusion  is  given  some  weight  by  the  fact  that  the  trend  curves  for  =1-5 

and  2. 67  show  that  the  effects  due  to  these  factors  are  small  at  y.  =  25.  Additional  data  points  at  the 
higher  mass  ratios  would  be  of  great  value  in  establishing  flutter  boundaries  at  higher  mass  ratios, 
however,  such  tests  were  beyond  the  scope  of  the  contract  and  were  not  feasible  with  available  models. 

c.  Data  at  OJ^/aJh,  =  2.67 

Data  obtained  at  this  frequency  ratio  are  shown  on  Figure  9.  The  flutter  boundary  for  /x  =  25 
is  defined  up  to  Mach  1.05  by  the  data  from  models  2  and  3.  A  sequence  of  run  points  obtained  for 
model  6  at  ^x  =  25  and  (0<£A)ti  =  2.64  shown  as  a  phantom  line,  has  been  plotted  between  Mach  1.05  and 
1.15.  During  this  sequence  the  model  was  observed  to  be  closest  to  flutter  at  Mach  1.09.  Examination 
of  the  oscillograph  records  indicated  that  the  model  had  not  quite  reached  flutter.  The  jj.  =  25  curve 
of  Figure  9  has,  therefore,  been  passed  above  this  test  sequence  line,  but  parallel  to  it  at  Mach  1.09. 
As  in  the  co^/tUh,  =  1.5  case  no  significant  transonic  dip  is  indicated.  Model  6  showed  a  decreasing 
flutter  tendency  above  Mach  1.09  indicating  a  definite  rise  in  the  fi  =  25  flutter  boundary  above  Mach 
1.09. 

Figure  9  shows  a  low  stability  region  as  a  dotted  area  between  Mach  1.0  and  1.2  belowthe 
^  s  25  flutter  boundary.  Mass  ratios  of  models  in  this  region  ranged  from  25  to  80.  The  jj.=  35 
point  for  model  10  at  Mach  1.18  was  obtained  by  holding  Mach  number  constant  and  slowly  increasing 
tunnel  density  (moving  up  on  Figure  9).  The  test  was  stopped  at  the  test  point  indicated  as  flutter 
seemed  imminent.  As  In  the  w^/Wh,  =1-5  case  raising  mass  ratio  above  25  thus  appears  to  drop 
the  flutter  boundary  (in  terms  of  the  parameter  Vg/br  oih,  )  at  low  supersonic  speeds. 

d.  Data  at  =  oo 

*  A  minimum  of  testing  was  done  with  the  pitch  degree  of  freedom  locked  out  in  view  of  the 
relatively  plentiful  amount  of  data  available  for  swept  cantilever  wings.  Consequently  only  two  data 
points  were  obtained  to  allow  some  correlation  of  pitch  flutter  data  with  this  limiting  case.  To  get 
flutter  in  the  desired  speed  range,  Mach  number  was  set  and  flutter  was  approached  by  raising  tunnel 
density.  Flutter  occurred  at  mass  ratios  of  34.8  and  35.7  at  Mach  numbers  of  1.0  and  1.06,  res¬ 
pectively.  The  points  are  plotted  on  Figure  10  in  terms  of  Vg/br  a>h|V^x.  This  figure  also  shows  the 
=  25  trendlines  for  ou,/cuh,  =  2.67  and  1.5  and  the  three  data  points  at  cu^/cof,,  =  0.9  for  comparative 
purposes.  The  drop  in  Ve/brCdf,,  accompanying  increasing  mass  ratio  above  25  which  has  been 
demonstrated  for  =  1.5  and  2.67  is  undoubtedly  the  reason  why  these  two  data  points  for 

0>(£/a;h|  =  ®  do  not  fall  above  the  cu^/cut,^  2-67  trend  line  for  /x  =  25. 

A  transonic  dip  followed  by  a  pronounced  supersonic  rise  in  Vg/brCOais  reported  in  Reference  3 
for  swept  cantilever  wings  of  aspect  ratio  3.1,  taper  ratio  0.4,  and  uin/u;a  =  0.2.  The  models  of  this 
program  do  net  seem  to  show  evidence  of  this  trend,  however  the  limited  cantilever  data  of  this  pro¬ 
gram  and  differences  between  models  of  the  two  programs  may  make  a  comparison  unjustified. 

3.  Effect  of  Pitch  Frequency  on  Flutter  Speed 

To  show  the  experimentally  observed  relationship  between  pitento  bending  frequency  ratio  with 
a  minimum  of  compressibility  influence,  all  data  points  up  to  Mach  0.92  were  plotted  in  Figure  13 
showing  Vf/bcCiJhi  versus  aiA/cdh,  ■  These  data  were  for  mass  ratios  between  23.5  and  24.6.  The 
drop  in  a)<£/a>h(  from2.67to0.9  causes  a  drop  in  Vg/br  cuh,  of  5%.  Analysis  values  of  Vr/brd>h|  versus 
from  Figure  6  are  also  shown  on  Figure  13  for  comparison.  The  analysis  curve  at  fJ.  =  25 
shows  a  drop  in  VR/br  0^,  of  6%  between  ou<£/Wh,  =  2.67.  and  0.9.  In  this  frequency  ratio  and 
Mach  number  range  the  experimental  reduced  velocities  average  7.5  percent  higher  than  analysis 
values. 

Figure  11  indicates  that  for  ou<f,/iOh|  between  1.5  and  2.67  the  effects  of  pitch  frequency  on 
Vg/br  0Uhi  are  not  significantly  alterea  by  compressibility  effects  in  the  transonic  range,  i.e.  the  jj.  - 
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25  curves  at  the  two  frequency  ratios  are  roughtly  straight  and  parallel  through  the  entire  transonic 
region.  From  Mach  0.95  to  1.20  the  drop  in  VE/br  uin,  in  going  from  o>h,  =  2.67  to  1.5  does  not 
exceed  7%.  Figure  6  shows  a  4.0  percent  drop  in  the  corresponding  "reference"  analysis  values,  in¬ 
dicating  a  general  correspondence  between  test  and  "reference"  analysis. 


The  effects  of  going  from  u>^/a>hi  =  2.87  to  ao  and  from  oj^/c o^=  1.5  to  0.9  can  not  be  accur¬ 
ately  established  due  to  the  simultaneous  variations  In  Mach  number  and  mass  ratio,  and  the  sparsity 
of  test  points  at  a»^/cohj  =  0.9  and  CO  .  Additional  testing  should  be  done  at  these  frequency  raHos 
holding  mass  ratio  constant  to  allow  quantitative  observations  of  the  separate  effects  of  pitch  to 
bending  frequency  ratio  and  mass  ratio  through  the  transonic  range. 

4.  Flutter  Frequency  Characteristics 

Flutter  frequencies  are  shown  ratioed  to  model  torsion  frequencies  and  analysis  flutter  fre¬ 
quencies  in  Figures  15  and  16  respectively.  Considering  the  ft  =  25  data,  a  drop  of  about  20%  occurs 
in  tuE  /aia  as  Mach  number  increases  from  0.75  to  1.2  with  the  greater  part  of  the  drop  occurring 
between  Mach  0.75  and  0.90.  This  would  indicate  a  more  prominent  participation  of  the  bending  mode 
as 'substantiated  to  a  limited  extent  by  the  amplitude  ratios.  Higher  mass  ratios  cause  a  progressively 
greater  drop  in  u;E/uJa  .  The  extreme  case  occurs  at  ai^/ajh(  =  0.9  where  the  combined  effects 
of  Mach  number  and  ft  increase  from  24.5  to  70.9  are  to  cause  a  drop  in  ajE/uja  of  45%. 

This  mass  ratio  trend  disappears  in  the  u)E/uiR  plot  of  Figure  16  indicating  a  reasonable  ana¬ 
lytical  prediction  of  mass  ratio  effects  on  frequency.  Analysis  appears  to  predict  the  proper  flutter 
frequency  in  he  vicinity  of  Mach  0.8  (<d£/a)R  =  1.0).  Since  the  analyses  performed  do  not  account  for 
,  compressibility  effects  (u>r  Is  not  a  function  of  Mr.ch  number)  the  drop  in  cuE/oia  observed  in  Figure 
.14  Is  repeated  in  Figure  15. 

5.  Amplitude  Ratios  and  Phase  Angles 

Both  analytical  and  experimental  amplitude  ratios  and  phase  angles  are  presented  in  Figure  17. 
As  may  be  seen,  analysis  predicts  a  phase  shift  of  <£,  a  and  ht  from  a  roughly  in  phase  condition 
relative  to  h,  to  a  nearly  out  of  phase  condition  as  drops  from  3.0  to  0.5.  Analytical  amplitude 

ratios  are  seen  to  increase  in  the  region  1 .0  i  and  then  to  decrease  at  very  low  oWcun,  with 

the  exception  of  4>/a  which  continues  to  increase  as  Oty/Wh,  decreases.  ' 

In  establishing  experimental  values  from  strain  gage  data  the  peaks  of  the  traces  were  examined 
first  of  all  to  ascertain  the  relative  phasing.  At  uM/a*,,  =  2.67  and  generally  at  uu/cuh,  =  1.5  the 
phasing  between  the  peaks  of  traces  was  either  negligible  or  was  of  the  same  order  of  magnitude  as  the 
estimated  reading  error.  Accordingly,  amplitude  ratios  were  established  for  most  data  points  assum¬ 
ing  zero  phasing.  For  model  4,  however,  at  10^/0^  =  1.5  and  for  the  three  data  points  at  co<£/tuh|  =  0.9 
phasing  between  peaks  appeared  to  be  significant  and  was  accounted  for  in  arriving  at  amplitude  ratios. 
The  general  Increases  in  a,  and  h£  relative  to  hj  as  ui^AUh,  drops  toward  0.9  indicates  a  replace¬ 
ment  of  the  bending  motion- by  a  more  prominent  pitch  motion  as  cu^/ojh,  drops. 
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Dotted  area  indicates  low  flutter  stability.  Phantom  line  is  fi  =  25  path  followed  in  getting 
flutter  with  model  6  at  Mach  1.17.  The  fi  =  25  flutter  boundary  must  pass  above  the  phantom 
line  for  lower  Mach  numbers  indicating  a  negligible  "transonic  dip”. 


CONFIDENTIAL 


Figure  B.Ve/brwh^  versus  Mach  number  for  W^Adh|  =  2.67  arid  p.  =  25  except  as  noted. 
Model  6  traced  out  the  phantom  line  at  ft  =  25  with  a  marginal  flutter  condition  at  Mach  1.09. 
The  M  =  25  flutter  boundary,  therefore,  passes  above  this  phantom  line  indicating  a  negligible 
"transonc  dip".  Dots  indicate  low  stability  area. 
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Figure  10.  Vg/brUif,,  -v^T  versus  Mach  number  showing  the  fj.  =  25  trend  lines  taken  from  Figures 
8  and  9  for  c oj.  /o)h  =  1.50  and  2.67.  Data  points  at  c /ajh  =  00  and  0.9  are  also  shown. 
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for  comparison. 


Figure  14.  Oscillograph  jccoru  of  flutter  of  model  6.  Test  point  115  of  first  test  period. 

=  1.53,  Mach  1.175,  fj.  -  28.0.  Model  diverged  in  flutter  about  ten  second 
after  record  was  taken. 
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Figure  15.  Experimental  flutter  to  torsion  frequency  ratio  versus  Mach  number, 
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SECTION  4 
CONCLUSIONS 


1.  For  a  mass  ratio  of  25  and  pitch  to  bending  frequency  ratios  of  1.5  and  2.67  the 

stabilizer  configuration  of  this  program  showed  a  negligible  sensitivity  to  compressibility  effects 
having  neither  a  pronouned  transonic  dip  nor  supersonic  rise.  For  oj^/o)h|  =2.67,  Vg/br  o;h|  V £T 
remained  between  1.74  and  1.80,  while  for  aj^/cuh|  =1.5  this  parameter  remained  between  1.64  and 
1.74  for  test  Mach  numbers  between  0.75  and  1.19.  For  all  of  the  above  test  conditions  the  ratio  of 
observed  flutter  speed  to  analytically  predicted  flutter  speed  remained  between  1.13  and  1.19. 

2.  The  conservatism  of  referer  je  analysis  appears  to  decrease  as  mass  ratio  is  increased  above 
25.  The  extreme  case  occureed  at  M  1.19,  cu^/u>h,  =  0.9  and  fji  =  70.9  where  the  analysis  gave  a 
flutter  velocity  10%  below  the  experimental  value.  It  appears  that  this  loss  of  conservatism  is  due  to 
the  relatively  small  Increase  in  flutter  speed  brought  about  experimentally  by  increasing  f-L  as 
compared  to  reference  analysis  predictions.  However,  further  data  are  required  to  definitely  est¬ 
ablish  that  compressibility  and  pitch  to  bending  frequency  ratio  are  not  significant  in  this  connection. 

3.  The  effects  of  variations  in  pitch  to  bending  frequency  ratio  with  mass  ratio  constant  at  25.0 
are  small  and  of  the  same  order  of  magnitude  as  are  predicted  by  reference  analysis.  Data  points 
below  Mach  0.92  at  /x  =  25,  show  a  drop  of  5%  in  Ve  /br  as  c o<f>  /ojf,,  drops  from  2.67 
to  0.9.  Analyses  show  a  corresponding  drop  of  8%.  The  drop  in\%/br  uj^  in  going  from  tu^/cuh|  =  2.67 
to  1.50  does  not  exceed  7%  at  any  point  between  Mach  0.8  and  1.15  at  fi  =  25.0. 

4.  A  region  of  low  stability  most  prominent  at  Mach  numbers  greater  than  1.0  exists  beneath  the 
flutter  boundary  at  all  pitch  to  bending  frequency  ratios. 

5.  Observed  flutter  frequencies  were  close  to  analyses  frequencies  in  the  vicinity  of  Mach  0.75. 
Increasing  Mach  number  from  0.75  to  1.2,  at  fj.  =25  resulted  in  a  drop  in  cog/oia  of  about  20%, 
indicating  Increased  prominence  of  the  first  bending  mode  in  the  flutter  mode. 

8.  Both  experiment  and  analysis  show  that  for  the  pitch  axis  location  used  the  pitch  degree  of 
freedom  does  not  become  significant  amplitude-wise  until  drops  down  close  to  unity. 

Analyses  show  that  the  first  bending  mode  shifts  from  a  roughly  in  phase  conditions  with  pitch  at 
ui^/ojht  =2.5  to  a  roughly  out  of  phase  condition  for  u)<p  /uih|  =1.0.  This  was  only  roughly 
confirmed  by  experiment. 
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SECTION  5 
RECOMMENDATIONS 


1.  The  models  of  this  program  were  designed  to  give  flutter  information  in  the  Mach  0.8  to  1.2 
range  at  fl  =  25  and  1.5  <  CU^/0Jh,<  4.0  •  At  the  time  of  testing,  interest  had  broadened  to 

include  down  to  0.9  plus  the  effects  of  increases  in  [i  .  The  very  limited  number 

of  flutter  points  which  could  be  obtained  per  model  prevented  complete  coverage  of  the  ranges  of 
variables  varying  only  one  parameter  at  a  time  to  allow  a  determination  of  the  separate  effects  of 
each.  It  is  therefore  recommended  that  additional  models  of  the  present  configuration  be  tested  at 
/J.  =25  for  cu^/u)h,  =  0.9  and  ®  for  the  Mach  0.8  to  1.2  range  to  complete  the  parameter 

study  at  ft  =  25.  A  repetition  of  much  of  the  /J.  =  25  program  at  /i.  =  45  and  70  would  be 
desirable  to  study  the  effects  of  mass  ratio  variations  particularly  at  -  0.9  where 

reference  analyses  apparently  become  unconservative. 

I.  In  view  of  the  Important  effect  of  pitch  axle  location  on  the  Influence  of  the  pitch  degree  of 
freedom  on  flutter  speed  which  has  been  demonstrated  by  WADC  tests  at  subsonic  speeds,  it  Is 
recommended  that  tests  be  performed  with  pitch  axis  as  a  variable  to  allow  the  effect  of  this  para¬ 
meter  to  be  established  at  transonic  speeds.  In  view  of  the  many  cases  of  all  movable  stabilizer 
flutter  which  have  occured  in  recent  years,  It  would  also  be  desirable  to  run  additional  tests  covering 
variations  in  aspect  ratio,  taper  ratio  and  sweep  to  allow  the  designer  to  avoid  critical  combinations 
of  flutter  parameters. 

3.  The  measurement  of  amplitude  ratios  and  phase  angles  during  flutter  Is  a  difficult,  costly  and 
time-consuming  operation  of  limited  accuracy  when  done  in  terms  of  Instrumentation  attached  to 
the  model.  It  Is  recommended  that  in  future  programs  high  speed  photographic  facilities  (2,000 
to  8,000  frames/ second)  be  developed  to  supplement  other  methods  in  the  study  of  flutter  modes  and 
general  model  performance.  At  the  present  time  camera  location,  loading,  triggering,  and  lighting, 
problems  usually  reduce  the  practicability  of  high  speed  photography. 

4.  The  present  program  illustrates  the  need  for  wind  tunnels  more  suited  to  flutter  testing. 
Desirable  features  of  such  a  tunnel  would  be:  (1)  blowdown  rather  than  continuous  flow  type  (or 
other  devices)  to  allow  higher  test  densities  (1-2  atmospheres)  to  avoid  the  danger  of  fan  blade 
damage  and  to  reduce  flow  turbulence;  (2)  many  access  doors  and  windows  to  allow  maximum  flexibility 
in  model  mounting  or  Injection,  model  clamping,  and  photgraphing  (up  to  8,000  frames/second)  with 
adequate  lighting. 
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APPENDIX  I 

Model  Details 


1.  Geometry  and  Dimensions 

A  sketch,  brief  description  of  the  test  configuration  and  summary  of  basic  characteristics  is 
given  in  Section  1.  Additional  dimensions  which  may  be  useful  for  reference  purposes  have  been  in¬ 
cluded  in  Figure  18.  Dimensions  shown  are  common  to  each  of  the  10  models  tested.  It  was  original¬ 
ly  Intended  that  the  elastic  axes  of  the  models  run  along  the  40%  chord  line.  Factors  of  unknown  origin, 
however,  caused  a  rearward  shift  resulting  in  an  average  EA  location  of  41.6%  chord. 

2.  Stabilizer  Construction 

The  models  were  of  single  beam  stressed  skin  construction  with  solid  styrofoam  core  used 
for  skin  stabilization.  Beam  width  and  skinthickness  dimensions  were  adjusted  from  model  to  model 
to  obtain  the  desired  range  of  frequencies.  The  construction  is  shown  in  Figures  19  and  20.  The 
main  beam  was  formed  from  7o  ST  aluminum  alloy  sheet  bent  into  two  channel  sections  and  riveted 
back  to  back.  Sheet  thickness  was  tapered  to  obtain  proper  mass  and  elastic  variations  down  the 
span.  The  I-beam  was  riveted  and  bonded  to  a  steel  root  fitting  between  stations  2.9  and  3.6.  The 
intersection  of  the  elastic  axis  with  the  fuselage  surface  is  designated  station  0.  To  provide  chord- 
wise  stiffness  0.020  inch  aluminum  ribs  were  spaced  at  2.9  inch  intervals  down  the  span.  For  Models 
1,  2,  5  and  6  which  were  tested  first,  flat  ribs  were  used.  For  the  remaining  models  the  ribs  were 
flanged  top  and  bottom  to  provide  greater  rib  stiffness  and  area  for  bonding  ribs  to  the  structural  skin. 

It  was  found  that  the  skin  buckling  loads  could  be  more  than  doubled  by  the  addition  of  light 
stringers.  These  were  formed  by  bending  0.004  inch  sheet  aluminum  into  angles  0.06  wide  and  0.10 
deep.  These  were  bonded  to  the  styrofoam  flush  with  the  surface  in  the  heavily  stressed  portions 
of  the  stabilizer.  Ballast  was  added  ahead  of  and  behind  the  main  beam  as  required  to  control  the 
c.g.  and  mass  moment  of  Inertia  of  each  model. 

The  skin  was  75ST6  sheet  ground  to  the  required  thickness  and  taper  to  provide  the  desired 
contribution  to  beading  and  torsion  stiffness  for  each  model.  The  sheet  was  bent  with  a  .016  radius 
(hot  formed)  around  the  leading  edge  and  bonded  closed  along  a  trailing  edge  seam.  The  skin  was 
bonded  to  the  internal  structure  using  carefully  controlled  quantities  of  bonding  resin  to  obtain  a 
satisfactory  bond  while  maintaining  accurate  control  over  mass  characteristics.  To  allow  transfer  of 
skin  loads  into  the  steel  root  fitting,  the  beam-root  fitting  joint  between  stations  2.9  and  3.6  was  con¬ 
toured  flush  with  the  akin  and  extended  in  the  chordwise  direction  as  indicated  in  Figure  19  to  provide 
a  large  bond  area  for  the  skin.  The  steel  root  fitting  mass  was  kept  closely  centered  about  the  pitch 
axis  to  avoid  excessive  pitching  moment  from  this  source. 

Bending  moment  sensitive  strain  gage  bridges  were  located  at  slut  Ions  1.2  and  11,5,  Tho  In¬ 
board  image  wan  muunied  la  (lie  eleel  *mmI  filling  while  Mm  millnmril  bridge  whs  IwmdPd  In  Hip  Inside 
of  the  stabilizer  skin  for  maximum  sensitivity.  A  torque  sensitive  strain  gage  bridge  was  placed  on 
the  steel  root  fitting  at  station  1.9. 

The  relatively  intricate  construction  used  for  these  models  was  considered  most  practical  to 
meet  the  jJ-  =  25  requirement  at  the  low  tunnel  densities  required  for  supersonic  speed,  and  at  the  same 
time  build  sufficient  strength  into  the  models  to  allow  them  to  resist  hours  of  exposure  to  low  sta¬ 
bility  and  high  turbulence  conditions  during  the  probing  for  flutter  boundaries. 
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(1)  E.A.  located  at  41. chord  measured  along  chord 
perpendicular  to  quarter  chord. 

2 

(2)  c.g.,  xa  and  ra  taken  in  direction  of  chord 

perpendicular  to  elastic  axis.  See  table  1.)  for  specific  values. 

(3)  MAC  based  on  span  to  center  of  fuselage  =  12.07" 

(4)  C/4  Taken  in  Windstream  Direction. 


Figure  18.  Stabilizer  Model  Dimensions  and  Geometry 
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BENDING  MOMENT  AND  TORQUE  SENSITIVE  STRAIN 
GAGE  BRIDGES 


BENDING  MOMENT 
SENSITIVE  STRAIN 
GAGE  BRIDGE 


LEAD  BALLAST 


-ALUMINUM 
I  BEAM 

SOLID  STYROFOAM  CORE 


ALUMINUM  RIBS 


Figure  19.  Typical  stabilizer  before  aluminum  skin  was  bonded  in  place.  Photograph  obtained 
with  rear  lighting  showing  internal  ballast  and  other  parts  agaii  st  translucent  styrofoar.,  core. 


Figure  20.  Crossection  of  Stabilizer  indicated  in  Figure  19 
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3.  Model  Mount  Characteristics 

Figure  21  shows  the  complete  test  assembly  mounted  in  a  jig  for  laboratory  tests.  The  luselage 
was  of  wood  with  the  center  section  hollowed  out  to  provide  space  for  model  control  mechanisms  as 
shown  in  Figure  22.  An  aluminum  cover  plate  was  used  to  close  this  section.  A  sponge  rubber  seal 
between  fuselage  and  stabilizer  wa3  provided  for  the  first  four  models  tested.  This  was  replaced  by 
a  close  clearance  contoured  block  (not  shown)  for  the  tests  of  the  remaining  models.  Pitch  excitation 
was  provided  by  an  electromagnetic  shaker  applying  an  oscillating  force  tangent  to  the  surface  of  the 
torque  tube.  The  lower  solenoid  shown  In  Figure  22  was  connected  to  a  linkage  which  locked  out  the 
pitch  degree  of  freedom.  The  upper  solenoid  was  connected  to  a  brakeshoe  type  damping  device. 

The  mechanism  for  controlling  pitch  frequency  is  shown  in  the  rear  view  of  the  assembly  ?.r> 
presented  in  Figure  23.  A  long  steel  spring  was  attached  to  the  end  of  the  stabilizer  torque  tube. 
The  carriage,  which  can  be  driven  the  length  of  the  spring,  clamps  the  spring  rigidly  by  means  of  a 
preloading  device  to  prevent  any  free  play  or  binding  due  to  variations  in  spring  thickness.  The 
carriage  track  rollers  are  preloaded  against  the  track  for  the  same  purpose.  Position  of  the  carriage 
was  indicated  through  a  revolution  counter  connected  to  the  carriage  drive  system. 

The  stabilizer  torque  tube  was  mounted  on  flexure  pivots  to  avoid  the  inertia  and  friction  of  bear¬ 
ings.  Strain  gages  attached  to  these  flexure  pivots  provided  pitch  displacement  information. 


Figure  21.  Test  assembly  mounted  in  laboratory  fixture.  The  circular  steel  plate  to  which  the 
fuselage  is  attached  is  installed  in  a  tunnel  wall  window  opening. 
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Figure  22.  Closeup  with  fuselage  cover  plate  removed.  Stabilizer  ro.-t  fitting  is  shown  bolted 
to  the  torque  tube  which  in  turn  is  supported  by  the  welded  hub  assembly  through  flexure  pivots 
The  torque  tube  extends  through  the  wall  and  is  attached  to  a  torsion  spring  as  shown  in  Figure 
23.  The  solenoids  at  the  left  operate  mechanisms  for  locking  out  the  pitch  degree  of  freetiom. 
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Figure  23.  View  of  pitch  restraint  system  behind  tunnel  wall.  The  pitch  restraint  spring  is  in  a 
vertical  piano  with  one  end  attached  to  the  stabilizer  torqu'  ube.  The  effective  length  of  the 
torque  tube  is  controlled  by  the  position  of  the  carriage  which  clamps  the  spring  firmly  while 
traveling  along  its  track  driven  by  the  two  screw  drive  shafts. 
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APPENDIX  n 

MODEL  CHECKS  AND  CALIBRATIONS 

1.  Model  Mass  Characteristics 

A  direct  measurement  of  the  distribution  of  weight,  c.g.  location  and  mass  moment  of  inertia 
could  not  be  made  by  slicing  the  models  up  after  wind  tunnel  testing  as  all  models  were  destroyed 
during  test.  However,  close  control  over  these  Items  was  maintained  by  measurement  of  the  weight 
of  each  component  prior  to  final  assembly.  The  proper  amount  of  ballast  at  each  station  was  then 
computed  and  Incorporated  into  the  models.  An  estimate  of  the  bonding  resin  used  to  bond  the  skin  in 
place  was  made  and  taken  into  account  in  computing  ballast  requirements.  Final  weight  and  Ia  dis¬ 
tributions  were  corrected  for  the  actual  amount  of  bonding  resin  used  by  assuming  that  the  difference 
between  actual  weight  and  design  weight  was  due  to  a  uniform  layer  of  bonding  resin  over  the  entire 
Inner  surface  of  the  skin.  (The  method  of  application  assured  this.)  This  procedure  left  a  minimum 
chance  for  error  and  provides  reasonable  assurance  that  the  design  values  of  c.g.  -  50%  chord  and 
r£  =  0.25  were  met.  Model  weight,  Iu  and  xa  over  the  span  are  given  for  each  model  in  Table 
m.  The  values  of  Xa  are  based  on  a  c.g.  at  50%  chord  and  elastic  axes  positions  determined  as  de¬ 
scribed  in  Appendix  n  2.  Mass  characteristics  for  the  non-uniform  station-widths  at  the  root  and 
tip  were  established  allowing  for  the  cutouts  due  to  sweep  in  those  regions.  Models  1,  2,  5  and  6  were 
ballasted  to  maintain  fi-  25.0  down  the  span  for  a  tunnel  test  section  density  of  0.001352  slugs/ft3. 
The  remaining  models  were  made  somewhat  lighter  by  using  the  minimum  ballast  for  each  model  that 
would  give  the  proper  non-dimensional  mass  characteristics. 

2.  Model  Elastic  Characteristics 

The  aspect  ratio  and  chordwise  stillness  of  the  models  were  sufficiently  high  to  allow  simple 
beam  theory  to  apply.  Model  El  and  GJ  distributions  and  elastic  axis  locations  were  determined  by 
load  deflection  measurements  for  each  model.  For  this  test  the  models  were  cantilevered.  Loads 
were  applied  at  station  17.4  on  the  40%chord  line,  and  10  Inches  forward  and  aft  of  this  point.  Verti¬ 
cal  translations  and  rotations  of  inboard  stations  were  measured  by  pairs  of  dial  indicators  as  shown 
in  Figure  24.  Dial  indicator  incremental  loading  effects  were  investigated  and  found  to  be  negligible. 
The  El  and  GJ  dictributions  of  Table  HI  were  found  by  adjusting  the  design  values  until  computed  de¬ 
flection  shapes  matched  the  experimental  data. 

Table  I  In  Section  1  gives  the  average  value  of  elastic  axis  location  determlnd  for  each  model. 
Since  the  c.g.  has  been  taken  at  50%  chord  for  all  models  the  elastic  axis  position  down  the  span  of 
each  model  may  also  be  seen  from  the  xa  distributions  given  In  Table  in. 

3.  Mass  and  Elastic  Characteristics  of  Pitch  System 

To  allow  any  desired  pitch  to  bending  frequency  ratio  to  be  realized  during  laboratory  and  wind 
tunnel  testing  it  was  necessary  to  determine  the  uncoupled  pitch  frequency  which  would  occur  at  any 
carriage  position.  Laboratory  calibrations  were  run  on  both  pitch  springs  for  this  purpose  using  a 
rigid  bar  mounted  to  the  torque  tube  in  place  of  a  stabilizer.  The  bar  was  given  a  pitching  moment  of 
inertia  of  21.27  lb.  in. 2  so  that  it  simulated  an  average  model.  Calibration  curves  of  pitch  frequency 
versus  carriage  position  were  made  for  each  spring,  and  correction  factors  were  established  for  each 
model  to  account  for  the  difference  in  pitching  Inertia  between  each  model  and  the  calibration  bar. 

The  effective  inertia  of  the  spring -torque  tube  combinations  were  determined  over  the  effective 
lengths  of  both  springs  by  vibrating  the  system  at  a  number  of  carriage  positions  with  and  without  the 
calibration  bar.  The  frequency  equations  were  then  solved  for  the  effective  mass  moments  of  inertia 
of  the  spring  torque  tube  combinations  at  the  respective  carriage  positions. 
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Figure  24.  Stabilizer  model  set-up  for  static  test  with  loading  bar  at  Station  17.4  and  dial 
indicators  at  four  inboard  stations  for  measurement  of  bending  and  twisting  displacements. 


Figure  25  shows  the  ratio  of  the  pitching  inertia  of  all  elements  of  the  pitch  system  behind  the 
fuselage  fairing  to  the  pitching  inertia  of  the  average  stabilizer  surface  outside  the  fuselage.  Figure 
25  contains  a  table  which  shows  the  factor  by  which  each  model  deviates  from  the  plotted  data  where 
minimum  weight  ballasting  was  employed.  For  the  pitch  to  bending  frequency  ratios  used  during  test, 
the  inertia  ratios  in  all  cases  remained  between  0.087  and  0.1171. 

To  determine  the  stiffness  of  the  torque  tube  support  system  against  roll,  a  rigid  aluminum  bar 
with  the  rolling  Inertia  of  the  stabilizers  was  cantilevered  from  the  torque  tube.  The  roll  frequency 
of  this  system  was  found  to  be  258  cps  with  an  effective  roll  axis  7.5  inches  behind  the  fuselage  fairing. 
Zero  air  speed  analyses  performed  with  and  without  the  roll  degree  of  freedom  showed  that  roll  had 
a  negligible  effect  on  the  stabilizer  model  coupled  frequencies,  but  did  influence  the  node  line  locations 
of  the  higher  mooes  of  vibration  to  some  extent.  Since  the  »-oU  frequency  was  roughly  2.5  times  the 
highest  flutter  frequencies  encountered,  roll  was  not  includ'  in  the  flutter  analyses. 
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Figure  25.  Ratio  of  the  total  pitching  inertia  of  all  elements  of  pitch  system  behind  fuselage 
fairing  to  the  pitching  inertia  of  the  stabilizer  surface  outside  of  the  fuselage  fairing 


4.  Model  Natural  Frequency  Characteristics 

a.  Computed  Uncoupled  Frequencies 

The  mass  and  elastic  characteristics  of  the  models  given  in  Table  in  were  established  by 
the  methods  indicated  in  Appendix  H— 1  and  II-2.  This  data  was  used  for  the  computation  of  the  un¬ 
coupled  bending  frequencies  using  Thompson's  method  and  uncoupled  torsion  frequencies  using  the 
Holzer  method.  The  uncoupled  frequencies  thus  obtained  are  listed  in  Table  I,  "Summary  of  Model 
Characteristics".  Bending  to  torsion  frequency  ratios  deviated  only  slightly  from  0.273. 

b.  Computed  Coupled  Frequencies  and  Node  Lines 

The  mass  and  elastic  data  established  in  the  foregoing  Sections  of  Appendix  n  were  em¬ 
ployed  in  zero  airspeed  analyses  over  the  range  of  pitch  frequencies  of  interest  to  establish  coupled 
frequency  ratios  and  node  lines.  The  first  and  second  bending,  torsion,  pitch,  and  roll  modes  were 
included  in  the  analysis.  Node  lines  for  the  first  three  modes  were  computed  at  Intervals  over  the 
range  of  pitch  to  bending  frequency  ratios  of  interest.  Analysis  node  lines  at  pitch  to  bending  fre¬ 
quency  ratios  of  0.9,  1.5  and  2.67  are  shown  in  Figure  26  as  dotted  lines.  Figures  27-29  show  the 
analytically  predicted  variations  in  the  coupled  ;.o  torsion  frequency  ratio  with  pitch  to  torsion  fre¬ 
quency  ratio  for  the  first  three  coupled  modes. 
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c.  Experimental  Coupled  Frequencies  and  Node  Lines 

In  the  laboratory,  node  lines  and  coupled  frequencies  for  the  first  three  modes  of  vibration 
were  obtained  over  the  range  of  available  pitch  frequencies.  The  observed  node  lines  for  Models  2, 
4,  6,  8  and  10  have  been  plotted  on  Figure  26  for  pitch  to  bending  frequency  ratios  of  0.9,  1.5  and 
2.67  as  solid  lines  for  comparison  with  the  analytically  predicted  values  shown  as  dash  lines. 

The  observed  coupled  frequencies  were  ratioed  with  the  uncoupled  torsion  fre  uencies  of  the 
respective  models  and  were  then  plotted  versus  pitch  to  torsion  frequency  ratio  in  Figures  27  -  29 
for  comparison  with  analytically  predicted  values.  The  general  agreement  between  analysis  and 
measured  characteristics  is  tr  i  n  as  substantiation  of  the  mass  and  elastic  characteristics  given  in 
Table  t.  I  and  HI. 

d.  Instrumentation  Calibrations 

Oscillograph  records  of  strain  gage  response  were  used  to  determine  model  natural  coupled 
frequencies,  flutter  frequencies  and  flutter  mode  shapes.  Bending  moment  and  torque  sensitive  bridges 
were  located  along  the  elastic  axis  of  the  model  and  a  strain  gage  bridge  was  placed  on  the  pitch 
flexure  pivot  springs.  The  latter  bridge  was  calibrated  directly  in  terms  of  pitch  displacement  of  the 
stabilizer  torque  tube.  The  bridges  on  the  models  were  calibrated  in  terms  of  bridge  response  to 
bending  moments  and  torques  at  the  respective  bridges.  Calibrations  were  also  run  to  determine 
bridge  temperature  sensitivity.  These  calibrations  allowed  the  determination  of  bending  moments 
and  torques  at  the  bridge  locations  on  the  models  during  flutter.  This  information  was  converted  to 
tip  amplitude  ratios  using  the  basic  methods  of  Section  A  3.13  of  WADC  TR  53-47  (Reference  2)  which 
assumes  that  the  flutter  mode  can  be  represented  by  the  superposition  (with  proper  phasing)  of  normal 

modes  of  vibration. 

« 
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Figure  27.  Ratio  of  coupled  to  first  torsion  frequency  vs  pitch  to  torsion  frequency  ratio. 
Solid  line  represents  analysis  results.  Dots  indicate  data  obtained  in  the  laboratory  for  each 
of  the  10  models. 
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APPENDIX  m 

FLUTTER  ANALYSES  DETAILS 


1.  Flutter  Equations 

Flutter  analyses  for  the  models  of  this  program  employed  the  methods  of  Reference  1  expanded 
to  four  degrees  of  freedom  by  the  addition  of  pitch  and  second  bending.  Using  the  notation  of  Figure 
28  the  displacements  h  (positive  down)  of  any  point  on  thq*  elastic  axis  and  rotation  a  about  the 
elastic  axis  (positive  with  leading  edge  up)  due  to  pitch  plus  the  elastic  degrees  of  freedom  will  be 


given  by  h  ^  <j>  y  sin  A-  K<£  +  fh|  h(  +  fhj,  h2 


and  a  r  cos  A  t  fag 

Substituting  these  values  of  h  and  a  for  the  corresponding  values  given  in  equations  (12)  of  Ref¬ 
erence  1  and  proceeding  as  usual,  the  flutter  determinant  elements  may  be  written  as  indicated  in 
Table  IV  where  the  flutter  equations  are: 


A  a  •>,  Ill  *  Aa*i+Aoh2!!l'0 

A*h,Ili  +  a*  +  A ^ ha hj.  =  o 

A.  .  hi  +  A,.  a  +  A.  .  <f>  +  A  h«-0 

h2hl —  hg  a  —  h2^  —  hghg_2, 
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The  aerodynamic  coefficients  Ach,  Aca,  Aoh  ,Aaa,  ACT  and  A0T  Included  in  the  determinant' elements 
are  defined  in  Reference  1.  Other  symbols  are  defined  in  the  list  of  symbols  at  the  front  of  the  re¬ 
port. 

2.  Mass  and  Elastic  Characteristics  Used  in  Analyses 

Flutter  analyses  were  run  in  dimensional  form  using  the  actual  mass  characteristics  of  Model 
1,  Table  in,  and  the  mode  shapes  shown  in  Figure  29.  The  desired  ranges  in  p.  and  Xa  were,  of 
course,  obtained  by  appropriate  variations  in  p  and  Sa  .  Since  the  models  were  constructed  and 
ballasted  to  have  common  non-dimensional  flutter  parameters,  the  characteristics  of  Model  1  were 
representative  of  the  whole  series.  Figure  25  shows  that  pitching  Inertias  of  those  elements  behind 
the  fuselage  fairing  varied  only  slightly  from  10  percent  of  the  stabilizer  pitching  inertia.  For  the 
analyses,  then,  the  total  pitching  inertia  was  taken  as  110  percent  of  the  average  stabilizer  pitching 
inertia  of  20.53  given  in  Table  I. 


Figure  29.  Mode  shapes  in  bending  and  torsion  used  in  flutter  analysis. 
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TABLE  IV 

DETERMINANT  ELEMENTS  USED  IN  FLUTTER  ANALYSIS 

V  '  Wf  [J.^Aekfw,  fh,dfh,/d]f  dj 

Ah|«2iS-fHftrd<J-Trf[l.b,Ac-fK)f^d^TANACb*ACTfKl  df</d}  dy.  ] 
A|,^*SfNAj[ln-fh|^^  +  COJAj’.Silfh1f*^  -K  X^X.^d?"  Ttf  [siNA^b'AeK^fh^Jf 
-SinAtamaX  ba(*dL)Acli  +  C05Aib1ACttfh|f^  d^-K4blAt[,fhif^  dy.] 

Ak.K,8  X*  fhffKidl^-¥p[Xb'A  chfh,fhx  <ty  -TAM  A  ft*  (*4J  Ack  fh)  d  fht  U)  dij  ] 

A«ht  1 X  S^ifh.d^-ITftX  b5A«Ki*fk,  dy  -TAN  a/.  b*(Vi*)  AaK  f«  d  fh,  /cty  d  ^  ] 

A**  8  (l  -  a)  /.  I  «f2  d.y.  -  TT  f>  [/.VA.«  fjd^  +TAN  A  X  b'  AcT  f*  df«  /dif  dy  ] 
A<+  1  SlNA/pSwf-f+y  dy  +  COSaX  Lf*fy  d^  -  K  Xif-^d^-TTf(siMAXl[b3A^ 

-  £(*/A  *) TANA  A«J f*f<f  dy  +  cosa X  A**  b^f*  dy-  K  X^3  A  aK  Xfy  d  ^  } 

4  SrfXifhjdy  “TT  p[X  b  A*k  f*  dlf  “ X b*(*^0  TAN  A  AakX  dXt  /d^  d^ 

A^k,  8  sinaXib  Vk,»  d^+  cos  A  X  S<  f4fh|  dy  -  K  X'm  ■ft  fh,  dy 

“FpfsiNAX  [  b*  Ack  fH(“  b*  TANA  Ack  dfh,/d^]  <f  d\f 

+COSA  rtbM  ck  ffc,  -  b4(‘/*„)TANA  A  ftk  cifk,  /dy  3  f*  d^. 

"K  X  [  bl  Ath  fK|-  b3  L*/&n)  TAN  A  Act,  ci fh,  /dij  ]  dt(.  } 

A*<8  SINA  X  d^  +  coSA/pI^f^f^  dif -KX^fyf*  dy 

“  TT  p  {  JIN  aX  tl>3Ac<f«t+  b^TAN  A  Act  df*/d|f  ]  }  f +  d  y. 

+ cosa(  [^AmX+^TAN  A  A*r  df«/d^dj-K£  [bM^f**  bVA«  A  ACTdC,Ajt^} 

A|4  8  [  I  -  A  ]  [  5INlA  £)n  f+'y  dy  +  l  SIN  A  COS  A  X  S*  if  ^  d  if 

+cos  A  X  I«  V  dy  -  2  K  sin  A  X  td  if  y  +  K*  X^  V  dif 

-2  K  COSaX  S.V  dif  ]-n'f{  SlNlA  X  [  blAck  y*-  bJ(^B)TANAAchJJ  ]f/d^ 

+3IN  A  COSa£  t (AaMc^-b^JlANA  A^Jf^dy+cosA  XVAm  if  d  if. 

-K SIN aX  [2^  b'Aeh-bs(*/4jTA«AAcJ f  * dtf-KcosAl  Aetf) ijdy\Cff%\>if d^} 
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1  SMA  /.hi  *  C03 A f#  S<f*f Kx  dy  -  K  fon  Vht  -  TTf  {sin  A  le  t  bl  Ath  f  K  v 

-  b*  (VaJtah  a  A«kdfhl  f*  a  y + COSA  jftb’A «nffcl-  b4  64)  tan  a  A  ^  dfhl  /J  *]f* 

-K/0  [blACKfKl- bl(*/iB)  TANA  Ach  ^ 

/\h, s  f^fs.dy -TTf  [/ob* Acs fh.fs,  d^  -TANA  fB bJ6/4  J  Aes  f  hl  dfK,  Aty  d  ^ 
Ahte<  8  /.^h^dy-ITf  t  /•*>*  Acif^f.,  dy  +  i  b4TAN  A  Act  f^dfrf/d*  dy 
Ahl*  s  SWAihifk^^  d^+C03A/XfKlf4  d  y. -Kj^Th fKtf* dy- rrp [smA {, b" A cs y  f h t dij. 
-SINATANACbY^^Acsfs^f^  d^4C03A/0b3AcAv^dj^  -  Kj,  b*Ats  fstff  dij.  ] 
[l-(wBt/fcj<)tn]4!hfK\  dy-TTf  [£ b'A  ch  fh‘  d^-TAN  A  /^(^JAtkf^Jfk/^ 
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